Asymmetric Wave Propagation Through Nonlinear PT-symmetric Oligomers 
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In the present paper, we consider nonlinear PT-symmetric dimers and trimers (more generally, 
oligomers) embedded within a linear Schrodinger lattice. We examine the stationary states of 
such chains in the form of plane waves, and analytically compute their reflection and transmission 
coefficients through the nonlinear PT symmetric oligomer, as well as the corresponding rectification 
factors which clearly illustrate the asymmetry between left and right propagation in such systems. 
We examine not only the existence but also the dynamical stability of the plane wave states and 
■ interestingly find them to be generically unstable. Lastly, we generalize our numerical considerations 

D ' to the more physically relevant case of Gaussian initial wavepackets and confirm that the asymmetry 

in the transmission properties persists in the case of such wavepackets, as well. 

o' 

INTRODUCTION 

p j Over the last fifteen years, and ever since its original proposal by Bender and co-workers the study of PT- 

symmetric Hamiltonian systems has become a focal point for numerous investigations at the interface between 
theoretical physics and applied mathematics. The fundamental appeal of such systems is that they respect key 
physical symmetries, namely the Parity (P) and Time-reversal (T) (but not necessarily the stronger constraint of 
the Hamiltonian being Hcrmitian), thus providing an intriguing alternative set of Hamiltonians with potentially real 
eigenvalues. In the context that we will examine below and for standard Schrodinger Hamiltonians with a complex 
potential V, the above constraints of PT symmetry amount to the potential satisfying the condition V(x) = V*(—x). 
£C) \ While this field commenced as, arguably, a mathematical curiosity associated with the foundations of quantum 
00 • mechanics, a major development from the point of view of applications arose recently. This was due to Christodoulides 
and co-workers [2| who realized that nonlinear optics (and the synthetic systems that can be engineered therein) may 
present a fertile playground for the experimental implementation of such PT-symmetric dynamics. Furthermore, the 
" added feature of nonlinearity typically present in such settings initiated the consideration of the effects of such PT- 
symmetric potentials on the nonlinear waves (such as bright or gap solitons) that may arise therein. Subsequently, the 
first realizations of PT-symmetry arose both in the nonlinear optics of waveguide couplers (i.e., either two waveguides 
. . . with and without loss |3( - the so-called passive PT- or in the more "standard" case of one waveguide with gain and 
one with loss Q) and also in that of electronic analogs thereof [f|. This progress and the perspectives that it enables 
towards future developments have, in turn, fueled a considerable volume of further theoretical studies both within 
the realm of linear PT-symmetric potentials 0- 14 1 and even in that of so-called nonlinear PT-symmetric potentials 



> 



(i.e., where the PT-symmetric prefactor appears in front of the nonlinear term) [l5- 17|. 



Another theme that has received considerable attention recently concerns the study of asymmetric (i.e., non- 
reciprocal) wave propagation in the context of various applications. Among the first examples discussed in the 
literature is the asymmetric phonon transmission through a nonlinear interface layer between two very dissimilar 
crystals (l8j . Later on, the theoretical proposition of a thermal diode (l9j induced relevant experimental realizations 
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In fact, recently this theme has 



ll| to act as unidirectional invisible 



in |2Cj ; similarly, an optical diode [21| was theoretically proposed 1221 a nd experimentally realized |23( . More recently 
similar proposals have been presented in left-handed metamaterials [2J] and have, in fact, been experimentally realized 
by different groups also in the context of acoustic waves in granular systems [2HI . 
been examined in the context of PT periodic structures which have been shown 

media (at least for sufficiently small extents of the periodic structure (27| ) with transmission coefficients and phases 
identical to the ones in the absence of the PT-structure. 

In the present work, we adopt a different perspective to that of these recent works and in fact one closer to the 
considerations of [28j| . The latter work considered a linear lattice in the presence of a set of N (embedded within 
the linear lattice) nonlinear elements. It was then shown that an asymmetry within the linear or nonlinear (still 
Hamiltonian though) properties of these elements would lead to an asymmetry of the relevant transmission of a plane 
wave through the chain. This was then extended to the more physically realistic case of a Gaussian wavepacket whose 
asymmetry was examined between incidence (on the asymmetric nonlinear region) from the left and from the right. 

Here, we consider these notions of asymmetric wave propagation but for PT-symmetric oligomers (namely dimers, 



trimers, quadrimers, etc.; see also the earlier study of such oligomers [12| which has motivated further recent studies 
such as the examination of quadrimers in [29^. In particular, we examine both the similarities and the differences 
with the above picture. In particular, we start our considerations by examining the transmission of plane waves 
through the chain. For such states, we analytically identify their reflection and transmission coefficients and obtain 
the corresponding transmittivity. We find some fundamental differences here in that such a quantity may exceed unity 
due to the presence of gain. Furthermore, we identify another fundamental difference in that beyond a particular gain 
strength a supercritical amount of "mass" may be collected at the gain site that may in turn lead to indefinite growth 
on this site. Nevertheless, for PT-coefficient strengths which yield meaningful transmittivities, a rectification factor 
is computed and clearly evidences the asymmetry between left and right propagation. Another aspect of the relevant 
plane wave solutions that we clarify concerns their generic dynamical instability, which we quantify and illustrate 
the (again, focusing) dynamical manifestation thereof. Lastly, we generalize our considerations (numerically) to the 
prototypical physically relevant case of an incident Gaussian wavepacket. We showcase that in that case as well, the 
transmission result differs between left and right such wavepackets impinging on the nonlinear PT-symmetric lattice 
region. 

Our presentation is structured as follows. In section II, we will give an overview of the relevant theory (some 
of the pertinent, more cumbersome details are relegated to an appendix). In section III, the analytical results are 
illustrated through numerical computations of both the transmittivities and rectification factors (both for the plane 
wave structures and for the Gaussian wavepackets) and of the existence, spectral stability and nonlinear dynamics 
of the plane waves. Finally, in section IV, we summarize our findings and present some directions for future studies. 



THEORETICAL ANALYSIS 



Stationary DNLS 

Starting our considerations with the existence of stationary solutions of the chain of interest, we examine the set 
of algebraic equations 

OJIpn = V n tp n - 4> n+ i - -071-1 + a n |l/> n | 2 l/>„ (1) 

on a one-dimensional lattice for ip n G C and oj € H. The parameters V n € C and are zero everywhere except 

for 1 < n < N so that the wave propagates freely outside of a finite region containing the nonlinearity and the 
(intended to be PT-symmetric, hence generally complex) linear potential. Rearranging ([!]) one obtains the so-called 
backward transfer map 

1pn-l = -1pn+l + (Vn ~ W + a n \lp n \ 2 )lp n (2) 

from which solutions to (TT]) can be constructed from knowledge of ip n at two nodes. 

The class of solutions whose (transmission and reflection) properties will be thercotically analyzed consists of plane 
waves of the form 

_ / R e ikon + Re~ ikon n < 1 , , 

Vn - | Telkon [A) 

with Rq,R,T G C representing the incident, reflected and transmitted amplitudes, respectively Note that outside 
of the nonlinear region, for n < 1 and n > N, ([3]) satisfies (fTJ) for any Ro,R and T if u = — 2cos(fco). Also directly 
from ([3]) [as applied to sites with n = and n = 1] , we have 

e- lk °iPo - 0! e' fc °0 o - 0i 
Ro = t r. — and R = —. -r- . (4) 

g — IKq _ glk glfe g — IKq v ' 

Thus, for any fixed values of ko and T, ipQ and ipi can be calculated by applying ([2]) iteratively starting with 
i(} N = Te ikoN and Vn+i = Te iko( - N+1 ^ from (0. Then, (gj gives the appropriate values of R and R so that ((TJ is 
satisfied at all nodes. 

For convenience we write the backward transfer map with n = N — I + 1 and in terms of *S> n for ip n d =' Te lkaN ^> n . 
This gives 

^n-i = + Sn-i+i^n-i+i (5) 

for Sj d =' Vj - uj + aj\T\ 2 \^>j\ 2 . For example, applying flSJ with N = 2 and * 3 = e ik " , * 2 = 1 gives 

S 2 = V 2 ~Lo + a 2 \T\ 2 (6) 

Si = Vi-Lu + ai \T\ 2 \8 2 -e lk °\ 2 
*o = -l + S 1 (S 2 -e iko ). 



Finally by (gj) and ([6]) we have 
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(7) 



and the corresponding transmission coefficient t d =' |T| 2 /|i?o| 2 is then 



t 



-,ik n 



-ik c 



l + (Si- e ife o)(e ife o _ S 2 ) 



(8) 



In the linear case (a± = a 2 = 0), it is immediately seen that t is the same for waves coming from the left or right 
side, independently of V n , as prescribed by the reciprocity theorem 28 . 3p| . 

It should be noted here that although the relevant calculation was presented for TV = 2, it can be performed 
for arbitrary values of N (naturally, the complexity of the intermediate steps is increased, the higher the value of 
N). For reasons of clarity of exposition but also of completeness, the relevant presentation has been given in the 
Appendix. We should also note that in comparison to the earlier work of [28|, there are fundamental similarities in 
the approach but also important differences in the nature of the results since our quantities Si (corresponding to the 
v and 8, respectively in 28|) are now, in principle, complex rather than purely real. Let us also indicate here that 
to quantify asymmetric propagation we will use the definition of a rectification factor / in the form: 



/ = 



t(kp,T) - t(-kp,T) 
t{k ,T)+t(-k ,Ty 



(9) 



whose non-vanishing values in the range 
metric transmission is tantamount to / = 



— 1,1] are measures of the asymmetry of transmission in the system (sym- 
0). 



Time Propagation 

We also briefly touch upon the tools that we will use towards the consideration of the dynamical evolution phe- 
nomena within our nonlinear Schrodingcr type chains 

i4> n {t) - V n <f> n (t) + 4> n+ i(t) + <f> n -i(t) = a n \(j) n {t)\ 2 (j) n {t). (10) 

In particular, in addition to direct numerical integration of Eq. (|10p . so as to monitor the dynamical evolution of the 
solutions, we will use a spectral stability analysis of stationary states (of the form ip n e~ lLJt discussed in the previous 
section) according to 

<f> n {t)= p n {t) + ep n {t). (11) 

Here, p n (t) = i\) n e~ %u)t and p n {t) = e~ lu)t (a n e wt + b n e~ w *) for w £ 1 and a n ,b n ,v G C. p n is assumed to be a 
standing wave solution of (TTU|) so that the resulting order-e equation for p n (t) is 

ipn - V n p n + p n+1 + p n -i + 0jp n = Q„ (2p n \p n | 2 + p\p* n ) . (12) 

The ensuing linear stability equations will yield the eigen-problem of the form: 



a n \ _ f Fi F 2 \ ( a r , 
K \ F 3 F 4 M 6* 



(13) 



for 



F 1 = diag(u - V n - 2a n \4' n \ 2 ) + G (14) 

F 2 = diag(-a n 4>l) 

F 3 = diag(a n (^* n ) 2 ) 

F A = diag{-uj + V* + 2a n \ip n \ 2 ) - G, 



where G is a sparse matrix with ones on the superdiagonal and the subdiagonal. Note that in (|1 3|) it is now 
convenient to think of a n and b n as column vectors. Given a stationary solution p n and values of V n ,a n which 
encode the nonlinearity for 1 < n < N, one then calculates the eigenvalues v in (|13|) . If v has a nonzero imaginary 
part this indicates that the perturbed solution (j) n (t) is unstable, as is easily seen by the form of p n {t) specified above. 
In practice, one diagonalizes a finite truncation of the matrix in (|13p . checking that the relevant eigenvalues are not 
affected by the truncation error. 
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FIG. 1: The figure presents the case of a PT-symmetric dimer with N — 2, ai,2 = 1, while Vi = iy and V% = —17. The left set 
of panels corresponds to the case of 7 = 0.25, while the right set of panels to the case of 7 = 0.5. Each set contains a contour 
plot of t(ko, T) (top left), a contour plot of f(ko,T) and a typical example of the dependence of t for ko — 2 (solid lines) and 
feo = —2 (dashed lines), so as to illustrate the asymmetry between the propagation for left and right incident waves. 
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FIG. 2: Same as the figure above but now for the trimer case of N = 3, 01,2,3 = 1, while Vi — ry, V2 = and V3 = —iy. The 
specific values of 7 for the panels given are 7 = 0.1 and 7 = 0.45. 

NUMERICAL COMPUTATIONS 

We now turn to computations in order to quantify the above theoretical results (as well as to extend beyond the 
range of what is analytically tractable). 

We start with the consideration of the transmittivity t(ko,T) and of the rectification factor f(ko,T) which are 
given for the case of the PT-symmetric dimer in Fig. Q] and for the PT-symmetric trimer in Fig. [2] The nonlinearity is 
uniform in both cases, but the linear potential is V\ = —V2 = i"/ in the former, while it is V\ = — V3 = 17 and V2 = 
in the latter. In both cases, a typical example of the transmittivity dependence for kg = 2 and kg = —2 is shown in 
the bottom panel of the figures. The asymmetry which is present in the top left panels between positive and negative 
values of kg and which is further quantified in the rectification factor of the top right panels clearly makes the case 
for the asymmetric wave propagation in these PT-symmetric oligomers. Although the values used here are below the 
PT transition of the underlying linear oligomer system, we have ensured that the relevant characteristic behavior 
(and presented asymmetries) exist both below and above that transition. However, as can also be inferred from 
the figures, the higher the value of PT-symmetric parameter 7, the stronger the manifestations of the asymmetric 
propagation of the waves. Another relevant observation to make here is that although the rectification factor is by 
construction bounded within [—1,1], the transmittivity is not bounded by unity (contrary to what is the case in the 
Hamiltonian example of [13 )• Hence, we can observe that in all the examples of t shown the relevant factor may 
exceed unity thus illustrating the existence of gain in the system. 

We now wish to touch upon one important aspect of the extended solutions ip n , identified in the form of Eq. ([3]) . 
Interestingly, a detailed examination of the stability properties of these solutions in fact reveals that they are gencri- 
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FIG. 3: The bottom panels show eigenvalues of the linearization of the example solutions shown in the top panel. In the 
PT-symmetric dimer (left), the parameters are N — 2 with ko = 2.5, T = 0.8, 7 = 0.1, while in the PT-symmetric trimer 
(right), they are: N = 3 with k = 1.1, T = 0.7, 7 = 0.1. 



cally unstable. A typical example of the corresponding spectral plane (Re(z/),Im(;/)) of eigenfrequencics v =Rc(^) +i 
Im(^) is shown for an example of the dimer and a corresponding one of the trimer in F ig. EH In fact, our computations 
indicate that this feature persists in the Hamiltonian variant of the model of Rcf. [28| (although a more detailed 
examination of the latter is of interest in its own right). The dynamical instability observed via the presence of 
imaginary eigenfrcquencies in the spectral plane of the PT-symmetric oligomers is corroborated by direct numerical 
simulations in Fig. HJ Here, we initialize the lattice with the solutions that were shown to be unstable via the spectral 
analysis of Fig. E3 to which a normally distributed small amplitude random perturbation has been superposed. The 
manifestation of the instability once again clearly underscores the role of the gain in the system. In particular, 
the dynamics reveals the tendency of the site which has gain to acquire a super-critical "mass" (or density), which 
subsequently grows exponentially well beyond the density of the linear background (or of the sites with loss). 

Given the above considerations and also the lesser physical relevance of exciting an initial condition of the form of 
Eq. (j3J) in an optical setting (which would be the prototypical potential realization of the considerations presented 
herein), we now turn to the examination of the dynamics of a more localized Gaussian wavepacket. The latter is 
assumed to be of the form: 

'iP n =Ie tkon - {n ~ no)2/s2 , (15) 

i.e., of amplitude /, centered at no G Z, while its speed is controlled by the parameter fco which assumes values 
in the interval — 7r < fco < 7r and s (chosen to assume the value 10 in what follows) is a parameter controlling the 
wavepacket width. It should be noted here that this aspect of our consideration is purely numerical (as it is not 
straightforward to obtain explicit analytical expressions for the transmission and reflection coefficients in this case). 
We thus evolve an incident wavepacket from the left or from the right and observe its propagation as illustrated in 
the space-time contour plots of Fig. [5] In what is shown below, we have performed the numerical computations for 
1=1, but we have confirmed that the relevant phenomenology persists for a wide range of non- vanishing values 
of I (approaching the linear limit of / — > approaches vanishing transmittivity differences, once again due to the 
principle of reciprocity). 

The transmission coefficient in this case is defined as the fraction of the total sum of the squared modulus of the 
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FIG. 4: The top panel shows the contour plot of the space n - time t evolution of \<j> n | for the solutions/parameters corresponding 
to the left and right, respectively, panels of Fig. [3] The rest of the panels show individual snapshots of the solution at t = 10, 40 
(for the left panels of the dimer case) and at t = 40, 80 (for the right panels of the trimer case). 



field across nodes that is transmitted to the right of the PT-symmetric oligomer after the Gaussian wavepacket passes 
over the relevant nodes. That is, we define t ko = J2 n <i \^n\ 2 /J2 n l^"l 2 for k o < and t ko = J2n>N l^nlVSn l^«l 2 
for fco > 0. We find that the transmittivity is greater for Gaussians approaching from the right and that the 
transmittivity difference is amplified as 7 increases, see Figure [S] which presents the relevant rectifying factor. This 
can be qualitatively understood once again on the basis of the structure of the PT oligomer. The node that has gain 
(which is on the left) favors reflection for a wavepacket from the left, while it favors transmission for a wavepacket 
impinging from the right, hence ensuring that i_fc > tk in our setup. Interestingly, in this case of the Gaussian 
wavepacket as well, for a high enough value of 7 the density accumulates and grows indefinitely at the nonlinear 
node that bears the gain. On the dimer with fco = n/2 this critical value is 7 ~ 0.6531, and with fco = — tt/2 it 
is 7 ~ 1.0890. On the trimer with fco = tt/2 the blowup occurs for 7 ~ 0.6737 and with fco = — ir/2 the value is 
7 ~ 0.8043. We note once again that these features are unique to the case of PT-symmetric oligomers through their 
gain-loss pattern and would be entirely absent in the earlier Hamiltonian installment of such chains in [28j . 
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FIG. 5: The space-time contour plot of |(/>n(i)| 2 is shown for left and right incidence on the PT-symmetric dimer (top panels) 
and on the PT-symmetric trimer (bottom panels). Top: dimer, N — 2 with no = —50, s — 10, fco = —tt/2 (left) and ko = n/2 
(right); Bottom: trimer, N = 3 with no = —50, s = 10, ko = — tt/2 (left) and ko = tt/2 (right), 
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FIG. 6: Rectifying factor for the incidence of a Gaussian wavepacket on the PT-symmetric dimer and trimer. Left: dimer, 
N = 2 with fc = tt/2, n = -50, s = 10; Right: trimer, N = 3 with k = tt/2, n = -50, s = 10. 



CONCLUSIONS 



In the present work, we considered a lattice setting where embedded in a linear Schrodinger chain was a nonlinear 
PT-symmetric oligomer, typically a dimer or a trimer. Our analytical considerations were focused around plane 
waves enabling us to analytically compute both the transmittivity and the rectification factor between left- and 
right-propagating such waves. These features amply evidenced the asymmetric nature of the propagation and even 
illustrated features particular to the gain-loss systems, such as the existence of over-unity transmittivitics. On 
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the other hand, we also considered the spectral stability of such states revealing their generic instability that was 
monitored and shown dynamically to lead to mass focusing on a single (gain) node of the lattice. This (as well as 
the more straightforward realization of a localized excitation in such lattices) , in turn, led us to the consideration of 
the asymmetry of propagation of a Gaussian wavepacket which we numerically quantified. Here, too, however the 
interesting phenomenon of the potential trapping of mass in a particular site with gain and the subsequent indefinite 
growth thereof were observed and quantified. 

These results suggest numerous interesting investigations for future work. It would be relevant to attempt to theo- 
retically quantify transmittivities and reflectivities of a Gaussian wavepacket perhaps through a judicious variational 
ansatz or some similar method that appropriately reduces the degrees of freedom, while taking into consideration 
both the complex PT-oligomer dynamics and the generic existence of a reflecting and a transmitting wavepacket. On 
the other hand, it would be particularly interesting to generalize the relevant considerations to higher dimensional 
settings and examine how incident waves from different directions may affect the transmittivity of different forms of 
two-dimensional PT-symmetric oligomers (in the simplest genuinely two dimensional case, PT-symmetric squares). 
These themes will be considered in future studies. 
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For the record we show five iterations of ([5]) starting from 'J/jv+i = e lk and ^jv = 1 



S N 

*N-2 

<5jV-2 
*JV-3 



= -*jv+i + S N ^ N 



(1 = 1) 



-e lk + 5 N 



V N ^ - u + a N ^\T\ 2 \S N ~ e lk \ 2 
-*JV + &V-l#i\T-l (1 = 2) 
-l + 8 N ^(S N -e lk ) 

V N - 2 -oj + ajv- 2 |T| 2 |l + £ w _i(e lfe - S N )\ 2 

-V N -l + <W-2*iV-2 (1 = 3) 

e ik -S N + 8 N _ 2 (-1 + 8 N ^(8 N - e* k )) 
Sn-2 + (e ik - 5 N )(1 - 6 N - 2 Sn-i) 

V N -3 ~ w + a7v-3|T| 2 |^w-2 + (<5iV - e ifc )(l - S N - 2 S N -i)f 

-*iV-2 + <5iV-3*iV-3 (i = 4) 

1 + ^-i(e Jfe - M + 5jv_ 3 (-5^-2 + (e lfe - <Sjv)(1 - S N - 2 Sn-i)) 
1 - <5jv-3<5jv-2 + (e* - 5n) (<W-i + 5 N - 3 (1 - S N - 2 5n-i)) 

Vn-4 —u + a7v-4|T| 2 |l - 8 N - 3 S N - 2 + (e* - 5 N ) (S N -i + 5 N - 3 (1 - S N - 2 5n-i)) | 2 

-^N-3 + S N -4^N-4 (1 = 5) 

= - (e* - 5 N )(l - 5 N - 2 S N -i) 

+5 N -4 (l - 8nsSn-2 + (e lk - $n) (Sn-i + 8ns(1 - 8 N - 2 8 N -i))) 
= 5n_ 2 + 6^-4,(1 — <5jv-3"5jv-2) 

+ (e lk - S N ) (Sn^^Sn^i + (8 N -4,8 N - 3 - - 8n- 2 8 N -i)) . 

Notice that if we iterate TV times we get S's equal to — u since V = a = for these nodes. Then Ro can be calculated 
from (j4]) in terms of appropriate 5's. 

For TV — 3, the algorithm gives a transmission coefficient 



#iV-3 



5n-4 

*JV-5 



g^™ g ^™ 



5i + (e i;fe -5 3 )(l-<52(<5i-e ifc )) 



(16) 



for 



S 3 = V3 - w + a 3 |T| 

5 2 = 7 2 -a; + a 2 |T| 2 |£ 3 -e <fe | 2 

£1 = yi-w + ai|T| 3 |l + $,(e <fc -*3)| 2 . 



(17) 



For JV = 4we obtain 



1 + <5 2 (e' lfc - 81) + (e lk - 5i)(5 3 + (1 - S 2 8 3 )(S 1 - e ik )) 



(18) 



for 



64 = V 4 -uj + a 4 \T\ 2 

S 3 = V 3 -oj + a 3 \T\ 2 \S 4 -e* k \ 2 

S 2 = V 2 - uj + a 2 \T\ 2 \l + S 3 (e lk ~ 8i)\ 2 

Si = Vi -uj + a 1 \T\ 2 \8 2 + (8 4 ~e tk )(l- 5 2 8 3 )\ 2 . 



(19) 



For TV = 5, we calculate 



g^™ g 



S 3 + (8! - e**)(l - 5 a * 3 ) + (e' fe - £ 5 )(5 4 (5i - e* fc ) + (1 - feWafa - e ifc ) - 1)) 



(20) 



for 
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S 5 = V 5 -uj + a 5 \T\ 2 (21) 

S 4 = V 4 -u + a 4 \T\ 2 \5 5 -e ik \ 2 

S 3 = V 3 - uj + a 3 \T\ 2 \l + 5^ - 5 5 )\ 2 

S 2 = V 2 - oj + a 2 \T\ 2 \5 3 + (6 5 - e lk )(l ~ 5 3 5 4 )\ 2 

Si = Vi - uj + ai\T\ 2 |l - 5 N - 3 8 N - 2 + (e lh - S N ) (<5jv-i + 5 N - 3 (1 - 5 N - 2 Sn-i))\ 2 ■ 



